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Abstract 

Galaxy surveys have shown that luminous galaxies 
are mainly distributed in large filaments and galaxy 
clusters. The remaining large volumes are virtually 
devoid of luminous galaxies. This is in concordance 
with the formation of the large-scale structure in Uni- 
verse as derived from cosmological simulations. How- 
ever, the numerical results indicate that cosmologi- 
cal voids are abundantly populated with dark matter 
haloes which may in principle host dwarf galaxies. 
Observational efforts have in contrast revealed, that 
voids are apparently devoid of dwarf galaxies. We 
investigate the formation of dwarf galaxies in voids 
by hydrodynamical cosmological simulations. Due 
to the cosmic ultra-violet background radiation low- 
mass haloes show generally are reduced baryon frac- 
tion. We determine the characteristic mass below 
which dwarf galaxies are baryon deficient. We show 
that the circular velocity below which the accretion 
of baryons is suppressed is approximately 40kms~^. 
The suppressed baryon accretion is caused by the 
photo-heating due to the UV-background. We set 
up a spherical halo model and show that the effective 
equation of state of the gas in the periphery of dwarf 
galaxies determines the characteristic mass. This im- 
plies that any process which heats the gas around 
dwarf galaxies increases the characteristic mass and 
thus reduces the number of observable dwarf galax- 
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1 Introduction: Formation of 
Structure in the Universe 

During the last couple of decades new extensive 
observations of the Universe were made using both 
ground-based telescopes and space instruments. 
These measurements have provided new insights into 
the structure of the Universe on various scales. A 
wide range of the electromagnetic spectrum emitted 
by cosmic objects has been studied. The wavelengths 
extend from very long radio wavelengths to energetic 
gamma rays. This observational progress has been 
accompanied by considerable effort in our theoretical 
understanding of the formation of different compo- 
nents of the observed structure of the Universe from 
small scales (galaxies and their satellites) up to the 
largest scales (clusters of galaxies and superclusters) . 

The standard picture of structure formation was 



suggested by White and Rees 70 more than 30 years 
ago. It suggests that gravitational instability drives 
the dark matter (DM) to cluster hierarchically in a 
bottom-up fashion. The gaseous baryons settle into 
to the DM haloes (namely, bound virial DM struc- 
tures) and via gas-dynamical dissipative processes 
cool, fragment and form stars. This simple picture 
prevails today, although our knowledge of the details 
and of the actual processes has evolved dramatically 
since then. The non-linear nature of the gravita- 
tional dynamics and the gas-astrophysical processes 
makes the problem of structure formation virtually 
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intractable analytically, and therefore the fiecp ..Id 
relies on numerical simulations. A substantial part 
of the progress in understanding structure formation 
in the universe is due to the increasing possibilities 
to make numerical experiments using the largest 
massive parallel supercomputers. In the eighties 
the best simulations handled 32^ particles whereas 
at present 1024'^ particles became a standard for 
numerical simulations, a increase of 2^^ roughly in 
agreement with Moore's law. 

The effort of observers and theorists brought about 
the so called standard cosmological model. This 
model is based on the idea that some kind of dark 
energy contributes about 70% of the total energy 
density of the spatially flat Universe in which the 
total energy density equals the critical one. The dark 
energy is responsible for the observed accelerated 
expansion of the universe. The simplest form of the 
dark energy is the cosmological constant, which was 
introduced already in 1917 by Albert Einstein in his 
paper about the cosmological solutions of the field 
equations of general relativity. The remaining about 
30% of energy density consists of matter. About 
85 % of this matter is made of unknown dark matter 
particles, the remaining 15 % is the contribution of 
"normal" baryonic particles well known to particle 
physicists. This means that the nature of more 
than 95 % of the matter in the Universe is not yet 
understood. 

According to this standard cosmological model, 
the main process responsible for the formation of 
the observed structures is gravitational instability. 
The initial seeds, which eventually became galaxies 
and superclusters and all other structures, came 
from the quantum fluctuations generated during the 
early inflationary phase: ~ 10~^^ sec or so from the 
beginning of the Big Bang. The power spectrum 
of these primordial fluctuations has been confirmed 
by measuring the temperature fluctuations of the 
cosmic microwave background radiation. These 
temperature fluctuations tell us the magnitude of 
the small fluctuations in the Universe about 300 000 
years after the Big Bang. 



One of the key features of the standard model is 
its simplicity. The expansion rate and the clustering 
properties are described by only few parameters 
which are measured at present with quite high 
accuracy. These parameters are the current rate of 
universal expansion, Hq = h x lOOkms"^ Mpc~^, 
the mass density parameter, ilmat, the value of the 
cosmological constant, the primordial baryon 

abundance, fib, and the overall normalisation of 
the power spectrum of initial density fluctuations, 
typically characterised by erg, the present-day r.m.s. 
mass fluctuations on spheres of radius 8h~^ Mpc. 

The initial power spectrum (created during the 
inflationary phase) does not contain preferred length 
scales. However, the horizon size at matter-radiation 
equality is mapped on this scale free spectrum. Since 
this scale (today around 100 Mpc) is well above the 
typical scales of observed objects structure formation 
is predicted to be essentially scale invariant: in a 
statistical sense the structures on scales of galaxy 
clusters are repeated on scales of galaxies. Typically, 
small objects merge together and form more and 
more massive objects. However, the small objects do 
not disappear within those larger objects but rather 
form a complex hierarchy of substructures. This 
hierarchical scenario predicts that our Milky Way 
galaxy is expected to have as many satellites (many 
hundreds) as a cluster of galaxies has galaxies. 
However, only a few dozen satellite of the Milky Way 
have been observed yet: this is the well known miss- 
ing satellite problem 33 43 . Other manifestations 



of the scale free power spectrum are the predicted 
large number of dwarfs in low density regions of 



the universe 17 49 or the predicted spectrum of 



mini- voids in the Local Universe [62| . 

A better understanding of the physics of structure 
formation on small scales, in particular of the correct 
modelling of baryonic physics, could solve these 
problems 22 . Based on semi-analytical models 



recently Maccio et al. [40] (see also Ref. 34 ) claimed 
that the long standing problem of missing satellites 
can be solved within the ACDM scenario. The 
basic idea behind this solution of the problem is 
that the haloes which host a galaxy of a given 
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measured rotational velocity are more massive than 
expected by the direct association of rotational 
velocity and the haloes maximum circular velocity. 
Since more massive haloes are less frequent the 
problem of missing satellites is solved. However, any 
non-baryonic physics that reduces power on small 
scales compared to the standard model will also 
improve the situation. It is well known that warm 
dark matter acts in this direction by erasing power 
at short scales due to free streaming 
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Altogether, we arrive at a picture in which dark 
matter particles form the backbone structure for 
all objects in the Universe from clusters of galaxies 
to dwarf galaxies. Normal matter (baryons) falls 
into the potential wells formed by the dark matter 
particles and forms the luminous objects. The details 
of this formation process must be followed using 
hydrodynamical simulations. Galaxy formation 
depends on many physical processes, some of which 
are very poorly understood, and some of which 
take place on sub-grid scales and therefore cannot 
be modelled fully and consistently and need to be 
fudged numerically by ad-hoc recipes. For example 
the formation of stars and the feedback of stars on 
the intergalactic medium take place on scales orders 
of magnitudes below the present day resolution. In 
cosmological simulations of galaxy formation typical 
masses of "stellar particles" are 10* h^^ Mq or more. 
Thus such a particle represents the formation and 
evolution of a large number of real stars. 

The paper is organised as follows: In Sec 2 we intro- 
duce briefly the method of cosmological simulations 
and discuss observational facts about voids and how 
they are defined in simulations. In Sec. 3 we present 
the baryon fraction in low mass dark matter haloes 
and introduce the characteristic mass. In Sec. 4 we 
introduce a spherical model for the halo gas and dis- 
cuss the relation between photo-heating of the halo 
gas and gas accretion. Our findings are summarised 
in Sec. 5. 



2 Cosmological simulations 

The cosmic microwave background (CMB) originates 
from the recombination of the hydrogen in the Uni- 
verse at redshift about 1000. Primordial small den- 
sity fluctuations are imprinted on the spectrum of 
the CMB temperature fluctuations which is closely 
related to the initial power spectrum of density fluc- 
tuations. The observed unique features of CMB tem- 
perature fluctuations constrain the cosmological pa- 
rameters as well as the normalisation of the power 
spectrum of density fluctuations. The most recent 
results of the WMAP experiment combined with the 
measurements of the baryonic oscillations and super- 



nova data 21 yield a Hubble parameter h = 7.01, 



the density of dark matter JIdm — 0.228, the bary- 
onic density fibar = 0.046, the cosmological constant 
r^A — 0.726, the normalisation of the power spectrum 
as = 0.82 and its slope n = 0.96. These parame- 
ters fix the cosmological model and therefore also the 
power spectrum of the initial perturbations. Starting 
with the initial conditions the simulations follow the 
growth of the perturbation in a universe that expands 
according to the chosen cosmological parameters. In 
this expanding universe the gravitational interaction 
is Newtonian. 

2.1 Initial conditions 

The first step of running cosmological simulations is 
to set up the initial conditions, namely amplitudes 
and phases of small perturbations at a very high 
redshift. Having in mind that the largest structures 
in the universe - superclusters and voids - have sizes 
of 10 to 50Mpc, the simulated volume should be 
significantly larger. However, we may be interested 
in the structure of a much smaller object such as 
our Milky Way galaxy or its satellites. In A^-body 
simulations each mass element is represented by a 
point-like particle. The mass resolution is limited by 
the total number of particles computers can handle 
at a given time. Thus, increasing the representative 
cosmological volume decreases the mass resolu- 
tion. To overcome this problem mass refinement 
techniques have been developed. Here we follow 



the algorithm proposed in Ref. 32 : To construct 
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suitable initial conditions, we first create a random 
realisation at the highest possible resolution. This 
depends on the available computers; at present we 
reach 4096'^ particles. The initial displacements and 
velocities of N particles are calculated using all waves 
ranging from the fundamental mode k = 27r/Lbox 
to the Nyquist frequency k^y = 27r/Lbox x 
Then the resolution can be decreased by replacing 
8* (i = 1,5) neighbouring particles by one particle 
with 8' higher mass which results in a distribution 
of 2048^, 10243, 5123, 2563 or 128^ particles. Using 
a smaller number of more massive particles we first 
run low-resolution simulations until the present 
epoch. In this simulation we select the regions of 
interest. Then we repeat the simulation but this 
time we preserve low mass particles inside the region 
of interest. Outside of this region we progressively 
replace small particles by massive ones creating 
shells of more and more massive particles until we 
reach the low resolution region of 128^ particles. 
This procedure ensures that our object evolves in 
the proper cosmological environment and with the 
right gravitational tidal fields. 

Based on the power spectrum of density fluctua- 
tions we set up the initial conditions of our simulation 
at an early redshift which depends on the size of the 
box and the required resolution. To move the parti- 
cles from the original Lagrangian point on a regular 
grid to their Eulerian position we use the Zeldovich 
approximation. The power spectrum of the generated 
Gaussian stochastic density field corresponds to the 
input power spectrum. The random nature of the re- 
alisation manifests itself on scales comparable to the 
box size where the spectrum is sampled only by a 
few modes. By this procedure the cosmological pa- 
rameters and the normalisation of the density power 
spectrum fully determine the initial conditions of the 
cosmological simulations, except of a random number 
which characterises the starting point of the random 
number generator. Thus, on scales comparable to the 
box size cosmic variance enters through this number 
which characterises a given simulation. On smaller 
scales - typically a quarter of the box size - the reali- 
sation follows very closely the input power spectrum. 
Starting from the initial matter distribution we sim- 



ulate the formation of cosmological structures. More 
precisely, we integrate the Poisson equation by the 
help of particle methods (Gadget and ART). 

2.2 The simulation method 

Our hydrodynamical simulations have been run 
with an updated version of the parallel Tree-SPH 
code Gadget 58 60 . The dark and baryonic matter 



distributions are represented by particles. The grav- 
itational forces are computed by a new algorithm 
based on the Tree-PM method which speeds up the 
force computation significantly compared with a 
pure tree algorithm. The hydrodynamical equations 
are solved by a Smoothed-Particle Hydrodynamics 
(SPH) method. The code uses an entropy-conserving 
formulation of SPH 58 which alleviates problems 



due to numerical over-cooling. 

The code also includes photo-ionisation and radia- 
tive cooling processes for an optically thin primor- 
dial mix of helium and hydrogen. For computing the 
thermal evolution it uses ionisation, recombination 
and cooling rates as given in Ref. |:29 . Since we will 
set up later in this paper an analytic model for the 
evolution of the intergalactic medium (IGM) we give 
here in detail the processes which affect the temper- 
ature, T, of the IGM. It primarily changes due to 
photo- heating, fi, and radiative cooling, F. More- 
over, it decreases due to the adiabatic Hubble ex- 
pansion, dT (X H dt, where H is the time-dependent 
Hubble constant. Cosmological structure formation 
changes the temperature adiabatically, dT oc dp. Fi- 
nally, changes of the chemical composition also affect 
the temperature, dT cx dS^X^, where Xi denotes the 
abundance of the atomic species i with respect to the 
cosmic baryon density, Xi — ni/n\,. In summary, the 
temperature evolution is given by 



dT = 



2 CH - F) At 



- 2HTdt 



3 fcenbSiXi 

2TdA Td(I],,Xi) 
^ 3 



(1) 



A 



YjrjXrj 



where A denotes the overdensity, A = p/{p). We 
take the atomic hydrogen and helium species, {HI, 
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HII, Hel, Hell, Helll}, into account. We assume 
a constant, primordial helium mass fraction of 



Yp = 0.24. 



The physics of star formation is treated in the code 
by means of a sub-resolution model in which the gas 
of the interstellar medium (ISM) is described as a 
multiphase medium of hot and cold gas 59 71 . Cold 



gas clouds are generated due to cooling and they are 
the material out of which stars can be formed in re- 
gions that are sufficiently dense. Supernova feedback 
heats the hot phase of the ISM and evaporates cold 
clouds, thereby establishing a self-regulation cycle for 
star formation. The heat input due to the super- 
novae also leads to a net pressurisation of the ISM, 
such that its effective equation of state becomes stiffer 
than isothermal. This stabilises the dense star form- 
ing gas in galaxies against further gravitational col- 
lapse, and allows converged numerical results for star 
formation even at moderate resolution. See [59j for a 
more detailed description of the star formation model 
implemented in the Gadget code. 

2.3 Dark matter haloes 

During the cosmological evolution dark matter 
haloes are formed which accrete more and more 
matter or merge with other ones. Within a simu- 
lation with one billion particles it is a challenge to 
find structures and substructures. To find structures 
at virial over-density one can use a friends-of-friends 
algorithm with a certain linking length (0.17 of the 
mean inter-particle distance for the ACDM model 
at redshift zero). The resulting particle clusters 
are in general tri-axial objects. Substructures can 
be identified as particle clusters at smaller linking 
lengths (higher over-densities). The more different 
linking lengths are used the better substructures will 
be resolved. Thus a whole hierarchy of friends-of- 
friends clusters have to be calculated. To this end 
we have developed a hierarchical friends-of-friends 
algorithm which is based on the calculation of 
the minimum spanning tree of the given particle 
distribution. We use also the bound-density-maxima 
(BDM) algorithm 31 1 which determines spherical 
haloes and their sub-haloes. The code removes 



unbound particles which is particularly important 
for sub-haloes. 

At all redshifts the sample of haloes is characterised 
by the mass function of the isolated haloes. This 
number density of haloes of a given mass depends on 
the power spectrum. A first very successful analyti- 
cal ansatz to predict the mass function of haloes has 
been made by Press and Schechter [52] . Later on this 
has been improved by Bond et al. pTand Sheth and 
Tormen [56] . The search for an accurate and univer- 
sal function which describes the number density of 
haloes found in simulation at different redshifts led 
to a set of fitting functions proposed by different au- 
thors |27l|63l[68l 



2.4 Voids 

With the first available large galaxy surveys it 
became clear that there exist large regions in the 
Universe which are not occupied by bright galax- 
Regions of all possible sizes devoid 
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les 

of galaxies can be seen in all redshift surveys. The 
observational discovery was soon followed by the 
theoretical understanding that voids constitute a 
natural outcome of structure formation via gravi- 
tational instability 23 48 . Together with clusters, 



filaments, and superclusters, giant voids constitute 
the large-scale structure of the Universe. 

More than 20 years ago, for the first time the 
sizes of voids have been estimated in different 
samples of galaxies with measured redshifts [7|. 
Later on voids in the CfA redshift catalogues were 



studied by 66 , in the Las Campanas redshift survey 
by [4^, in the Optical Redshift Survey and in the 
IRAS 1.2-Jy survey by (9)[l0] and in the PSCz 



catalogue by Plionis and Basilakos 50 and Hoyle 



and Vogeley 25 . For a review of early observational 
efforts see Peebles 49 . More recently voids have 



been studied by many authors using the Two-Degree 
Field Galaxy Redshift Survey (4l|6l[26l|47l and the 



Sloan Digital Sky Survey 16 24 
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There were several attempts to find dwarf galaxies 

The overall 



in few individual voids 19 35 38 51 
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conclusion is that faint galaxies do not show a strong 
tendency to fill up voids defined by bright galaxies. 
The strongest arguments that voids are not popu- 
lated by dwarf galaxies were given by Peebles 49 



who points out that the dwarf galaxies in the ORS 
catalogue follow remarkably close the distribution of 
bright galaxies: there are no indications that they 
fill voids in the distribution of bright galaxies. To 
summarise, observations indicate that large voids 
found in the distribution of bright (~ M*) galaxies 
are empty of galaxies, which are two magnitudes 
below M*. Recently, Tikhonov and Klypin [62] 
studied the properties of mini-voids in the nearby 
universe. 



The void phenomenon was already the target of 
many theoretical studies [l]|8 12-15 36 
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Using a set of DM only simulations 17 have studied 
the inner structure of voids. The haloes in voids are 
arranged in a pattern, which looks like a miniature 
Universe with the same structural elements as the 
large-scale structure of the galactic distribution of 
the Universe. There are filaments and voids; larger 
haloes are at the intersections of filaments. The 
mass function of haloes in voids is much steeper 
for high masses resulting in very few galaxies with 
circular velocities Ucirc ~ lOOkms^^. Note, however, 
that in DM simulations it has been usually assumed 
that each dark matter halo hosts at least one galaxy. 
This may not be true. Physical processes of galaxy 
formation are not well known and there could be 
processes that strongly suppress the formation of 
stars inside small haloes, which collapse relatively 
late in voids. This issue will be discussed in detail in 
the following sections. 

There is a lot of interest in the void phenomena. 
However, there is no clear definition of a void. 
Almost every study uses its own definition of a void 
and its own void finder. Recently, the different 
approaches have been compared within the Aspen- 
Amsterdam Void Finder Comparison Project [5]. 
Here we use a rather simple and direct definition. 
In order to identify voids, we start with a selection 
of point like objects in 3D. These objects can be 
haloes above a certain mass or a certain circular 



velocity or galaxies above a certain luminosity. Thus 
the voids are characterised by the threshold mass 
or luminosity. The void finding algorithm can be 
applied both to numerical and observational data. In 
numerical simulations it takes into account periodic 
boundary conditions. At first it searches the largest 
empty sphere which is completely inside the given 
volume. The radius of this empty region is i?void- 
Then we extend this region to an empty region with 
irregular shape. We define this empty region as 
union of all spheres with radius Text > 0.789 i?void 
(thus with a volume greater or equal half of the 
volume of the original void sphere) which can be 
moved to their place without crossing any object or 
the box boundary. To find the other voids we repeat 
this procedure however taking into account the 
previously found voids. We now search the largest 
empty sphere which is completely inside the volume 
and does not intersect with any already known void. 
The voids are characterised by volume, shape and 
orientation. 

For the re-simulation we selected a void region from 
a larger periodic computational box. To construct 
suitable initial conditions, at first we replaced in a 
box of side-length 50 h"^ Mpc with 2048^ particles 
the low mass particles by 256'^ massive ones. After 
running the simulation at this low mass resolution 
we determined a void region of interest and, finally, 
we were rerunning the simulation with the original 
set of particles inside the void (2048"^) and replacing 
only outside the void the low mass particles by more 
massive ones. 



The baryon content of dwarf 
galaxy haloes 



The number density of dwarf galaxies observed in 
the Universe is apparently smaller than the predicted 
number density of low mass dark matter haloes 
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The missing galaxies are either not formed because 
the corresponding dark matter haloes do not exist, 
or no stars are formed in small haloes. As shown by 
Hoeft et al. f22^ and Okamoto et al. [46^, the op- 
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refined 
region 


simulation 
bos size 


mass resolution 
(dark matter) 


void 
filament 


50 ft-^Mpc 
80 h-^Mpc 


1.03 X 10^ Mrs, 
8.24 X 10^ h'^ Mq 



Table 1: Main characteristics of the two simulations 
analysed for determining the characteristic mass. 



tically thin UV-background already causes baryon 
deficiency in small haloes. In this section we de- 
termine the mass scale below which photo-heating 
causes baryon deficiency. 

3.1 Mass accretion history 

The gravitationally bound structures, i. e. galaxies 
and clusters of galaxies, grow in two ways: they 
accrete surrounding matter and they merge with 
other galaxies or clusters. Numerical simulations 
allow us to follow the merging history of the haloes. 
Every massive halo at the end of the simulation has 
a large number of progenitors which subsequently 
merge to the final structure. 

We compute the mass accretion history for all 
haloes in our two simulations, see Tab. 1. To this end 
we determine the haloes in each simulation snapshot. 
Then we go backward from the final snapshot and 
search for the most massive progenitor in the previ- 
ous snapshot. Since we have stored more than one 
hundred snapshots for each simulation we can simply 
search for the most massive progenitor in the vicinity 
of a given halo. However, during merger events the 
halo finder sometimes fails to identify the halo at all, 
since two approaching haloes appear spuriously grav- 
itationally unbound. Fig. [l] shows the mass accretion 
histories of ten sample haloes, where the jumps are 
caused by merger events. 

3.2 Characteristic mass 

We determine for all distinct haloes in our simulations 
the baryon fraction within the virial radius, /b = 

Mbar/Mtot = (Metar + Mgas)/(Mstar + Mge^ + Mdm), 

where Mbar denotes the baryonic mass, Aftot the to- 



tal mass, Mdm the dark matter mass, Mgas the gas 
mass and Mgtar in stellar mass. We take those haloes 
into account which consist of more than 180 dark 
matter particles. The resulting lowest halo mass is 
1.9 X 10^ h-^MQ and 1.5 x 10^ h'^ Mq for the 
void and for the filament region, respectively. More- 
over, we exclude all haloes which have been within 
the virial radius of a more massive halo in the past, 
because the dark matter and baryonic masses may 
be affected by stripping. Fig. [2] shows the resulting 
baryon fraction for both the low density region sim- 
ulation the filament region simulation. The low mass 
haloes evidently have a baryon fraction much smaller 
than the cosmic baryon fraction. Since the transition 
from 'cosmic mean baryon fraction' to 'baryon defi- 
cient' occurs within a rather small mass range this 
phenomenon can be characterised by a single param- 
eter, namely the mass at which the baryon fraction 
amounts in average to half of the cosmic mean, Mi/2- 
We determine the characteristic mass computing the 
average baryon fraction in mass bins and approxi- 
mating the resulting curve by the following analytic 
expression 



/b =(/b) 



1 ^ /logMtot 
— arctan 



log Af 1/2 



(2) 



22 



to 



We use a different expression than in Ref. 
make the fitting procedure more robust, in particular 
when a second free parameter is included, here the 
width wm- Moreover, we expect that the halo gas, 
which is present independently from the halo size, 
leads to a fioor in the baryon fraction. As we have 
shown in Ref. 



22 



the characteristic mass depends 
very little on the mass resolution used in our sim- 
ulations. This indicates that the physical processes 
causing the baryon deficiency are reasonably resolved 
even if there are only a few tens of gas particles in 
small mass haloes. 



In Ref. 22 



we focussed on the baryon content 
of dwarf galaxies in cosmological voids. The result 
we obtained might be affected by restricting to 
void dwarf galaxies. To investigate also the effect 
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Figure 1: The mass accretion histories for several haloes. In each panel the total mass (solid line) and the condensed 
mass (dashed line) of one halo is shown. We have divided the condensed mass by the mean cosmic baryon fraction, 
(/b), to highlight when the condensed mass does not follow the expectation according to the mean baryon density. 
In addition, the characteristic mass (dotted line) is shown as derived from baryon fraction in newly accreted matter. 
One can clearly see that haloes cease to condense baryons when their total mass falls below the characteristic mass. 



of the cosmological environment we include here 
dwarf galaxies in a cosmological filament. These 
galaxies have a higher merger rate than void dwarf 
galaxies with a similar mass but they do accrete 
less diffuse gas, see Ref. 11 and references therein. 
We have excluded all galaxies which have passed 
through the halo of a more massive galaxy since 
their dark matter and baryonic mass could be 
affected by stripping. Despite the cosmologically 
quite different environments, the characteristic mass 
of both regions coincide well, see Fig. [2] (top row). 
In conclusion, the general mass accretion history 
of a dwarf galaxy has only little effect on their 
baryon content, i. e. there is universal characteristic 



mass for all dwarf galaxies. Note, this result is 
obtained with a homogeneous UV-background. We 
will argue later that the UV-heating may have a cru- 
cial impact on the baryon accretion in dwarf galaxies. 



Observationally, the total mass of a galaxy is 
difficult to determine. Galaxies are preferably 
characterised by the maximum circular velocity, 
Wmax- One should note that also i^max is difficult to 
measure in dwarf galaxies since most rotation curves 
do not approach a maximum, i. e. the maximum of 
the circular velocity lies outside the observable gas 
distribution. However, we will characterise haloes 
mainly by Wmax- Fig. [5] (middle row) shows the 
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Figure 2: The baryonic mass fraction in units of the cosmic mean, Qh/^m- Filled squares refer to the simulation 
of low density region while stars refer to the filament region. Only haloes with more than 180 dark matter are taken 
into account. The top row shows the baryon fraction at different redshifts as a function of total halo mass. The 
middle row gives the baryon fraction as a function of maximum circular velocity. The bottom row shows the baryon 
fraction of newly accreted matter (see text for a more detailed description). 



baryon fractions in distinct haloes as a function of 
the maximum circular velocity. 

Our aim is to determine what makes haloes below 
the characteristic mass baryon deficient. The baryon 
content at a given time is accumulated over the his- 
tory of the halo. To assess what prevents at a given 
time haloes from accreting baryonic mass we deter- 
mine the baryon fraction of newly accreted matter 

j^accr ^ Mcond{ti+l) " ^4ond(^i-l) 

" Mtot(i,+i)-Mtot(t.-i) ' ^ ^ 



where and ti^i denote the masses of a given halo 
at an earlier and at a later time, respectively. Fig. [2] 
(bottom row) shows that small haloes do not accrete 
baryonic mass. Since we consider the differential evo- 
lution of the dark matter and the baryonic masses 
we have a large scatter in the characteristic mass 
of newly accreted matter. However, similar to the 
baryon fraction in haloes there is a clear separation 
between low mass haloes which accrete virtually no 
baryonic matter and high mass haloes which accrete 
matter with the mean cosmic baryon fraction. 
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Figure 3: Left •panel: Evolution of the characteristic mass of the baryon fraction in distinct haloes (squares) and in 
newly accreted matter (circles). For comparison the mass accretion histories of five haloes, as depicted in Fig.[lJ are 
shown. Right panel: Evolution of the characteristic circular velocity in distinct haloes and in newly accreted matter. 
Additionally, the circular velocity histories of the haloes in the left panel is shown. 



3.3 Evolution of the characteristic 
mass 

Fig. [2] shows that the characteristic mass decreases 
with redshift. For instance, at redshift z ~ A the 
characteristic mass is about 3 x 10* Mq, while 
it amounts to 6 x 10^ Mq at ^ = 0. Fig.|3] 
shows the redshift evolution of the characteristic 
mass as obtained from combined halo lists from our 
two simulations. For comparison also a few mass 
accretion histories from Fig. [T] are shown. 

Plotted as a function of circular velocity the 
evolution of the characteristic mass of the newly 
accreted matter has a rather simple behaviour: after 
a steep increase it is almost constant at ~ 40kms~^ 
for z < 3. Note, in Ref. 22 we found a lower 



value since we used the circular velocity at the 
virial radius, i'(rvir), instead of Wmax- Moreover, 
we considered only the baryon fraction in haloes. 
Since semianalytic models of galaxy formation 



assume a low-mass cut-off for gas accretion, see 
e. g. Ref. 57 and Ref. 37 , we investigate here 
also the characteristic mass of newly accreted matter. 

Because the circular velocity of most of the dwarf 
galaxy haloes is almost constant with time, except a 
major merger occurs, one can separate two types of 
evolutionary scenarios for dwarf galaxies: (i) small 
galaxies with v^ax ^ 40kms~^ accrete baryons only 
for 2; > 4 and (ii) galaxies with v^a^ ^ 40kms~^ are 
not affected at all. In Fig. [l] we have also depicted for 
each halo the evolution of M^J^'^ as shown in Fig.jSj 
(left panel) . One can see nicely that haloes get baryon 
deficient when their halo mass is below My" . 

3.4 Can the characteristic mass ex- 
plain the void phenomenon? 

The void phenomenon comprises two observational 



findings 49 : First, the distribution of massive galax- 
ies (> i*) shows large empty regions, the 'voids', see 
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Sec. 2.4. Low-luminosity galaxies populate the voids 
but with a low volume density compared to the dense 
Secondly, the galaxies found in voids 



24 



regions 

are not special at all. According to the concordance 
ACDM cosmology voids are filled by dark matter 
haloes with the size of dwarf galaxies 17 . To judge 



if this is at odds with the void phenomenon in ob- 
servations one has to assign a luminosity to the dark 
matter haloes. Recently, Tinker and Conroy [64] 
applied a Halo Occupations Distribution (HOD) to 
a large dark matter simulation and analysed void 
regions. They found that the haloes in void regions 
which host galaxies with an absolute magnitude 
Mr > —10 are located close to the walls of the 
voids. Only even fainter galaxies, i. e. galaxies 
which are significantly below the current detection 
limit of surveys covering large void regions 24 , 
populate the entire void volume. According to 
Tinker and Conroy galaxies with Mr — —10 reside 
in haloes with a mass ~ 10^^ h~'^ Mq. Hence, 
the photo-heating discussed here would reduce the 
luminosity of the haloes filling the entire void volume 
rendering those galaxies even more difficult to detect. 

The HOD used by Tinker and Conroy implies that 
the mass-to-light ratio increases by a factor of ten 
from haloes with a mass of 10^^ Mq to haloes 
with 10^^ Mq, independent of the environment. 
This is consistent with the compilation of the baryon 
content of galaxies by McGaugh 42 . Haloes with 



the mass of 10^° Mq seem to host only 5 % of the 
baryonic mass corresponding to the cosmic mean. 
In our simulations we do not find any indication 
that the UV-background reduces the baryon fraction 
in galaxies significantly more massive than the 
characteristic mass. Therefore, we speculate that 
stellar feedback plays a key role in reducing the 
baryon fraction galaxies with Mtot > Mii2- Note 



that the baryon fractions in Ref. 42 shows indeed 
a break at Vd^c ~ 30kms~^, consistent with our 
findings for the characteristic mass. 

In summary, the surprising low number of low- 
luminosity galaxies in voids can be explained by a 
mass-dependent mass-to-light ratio, as obtained from 
a global Halo Occupation Distribution and from a 



proper mass modelling of galaxies. Haloes above 
the characteristic mass are expected to reside close 
to the walls of voids. The photo-heating discussed 
here reduces the luminosity of the low-mass and low- 
luminosity galaxies which can populate the entire 
void volume. Current surveys are limited to detect 
galaxies close to the walls due to the sensitivity limits 
of the surveys. 

4 Why small haloes fail to ac- 
crete gas 

In the previous section we have shown that dark mat- 
ter haloes with Wmax ^ 40kms~"^ do accrete dark 
matter for z < 4 while they do not accrete baryons for 
smaller redshifts. Halos with Wmax < 35kms~'^ are 
generally baryon deficient for z — {). These results are 
obtained for high-resolution hydrodynamical cosmo- 
logical simulations including a homogeneous photo- 
heating of an optically thin medium and including 
radiative cooling for a primordial chemical composi- 
tion. In this section we investigate in detail the origin 
of the baryon deficiency which is evidently related to 
the photo-heating by the UV background. To this 
end we set up an analytic model and investigate how 
baryons cool in dwarf galaxies. 

4.1 Outline of the analytic model 

Our simple model is constructed as follows: In the 
beginning the gas in the universe is homogeneously 
distributed. In the process of structure formation gas 
is accreted by galaxies, i. e. its density increases. 

Starting from a homogenous distribution the gas 
in the periphery of haloes and in the IGM follows 
the evolution of the dark matter. For our model 
here we assume that there is universal evolution of 
the overdensity, A(t, Aq), which only depends on 
the time and on the final overdensity, see Fig. [4] By 
the help of Eq. we can determine the tempera- 
ture as a function of the overdensity. For a given 
temperature-density relation we can determine the 
gas density profile of NFW-type haloes. Assuming 
that the gas density in the periphery of the halo is 
given by the mean cosmic baryon fraction we can 
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Figure 4: The overdensity evolution A(a,Ao). We 
group in the simulation snapshot for 2 = all particles 
according to their overdensity. Then trace back the den- 
sity of the SPH-particles and compute the average den- 
sity of the groups for earlier redshifts (solid lines) . More- 
over, we show three approximated overdensity evolutions 
(dashed lines) according to Eq. Q. See Sec. 4.2 for de- 
The peak-like structures at z ~ 6 are caused by 



tails 

photo-evaporation and are neglected in our approxima- 
tion. 



determine the central gas density and temperature. 
If the resulting cooling time is short the halo 
would accrete baryonic matter, if the cooling time 
is long the halo would fail to accrete baryonic matter. 



4.2 Effective equation of state of the 
IGM 

For a known overdensity evolution, A(i, Aq), we 
can integrate the thermal evolution of the IGM, 
Eq. ([1]). To determine the heating and the cooling. 
Hit) and T{t), respectively, we have to know the 
abundance of each species, i. e. the ionisation state 
of hydrogen and helium. As in the cosmological 



simulations we adopt collisional equilibrium, hence, 
we can determine for a given time and a given 
temperature the relative abundances. We solve 
the rate equations for collisional ionisation, for 
recombination, and for photo-ionisation. Now we 
can compute the photo-heating and cooling of the 
gas. Time-dependent specific energy injection rates, 
Ei, are tabulated according to the UV-background 
model. As a result the total photo-heating can 
be written as H = nnieHi + ^T-HeieHei + "HeiiCHeii- 
Cooling rates are used as given in [29]. Hence, 
for a given overdensity evolution, A(t, Aq) and 
UV-background model we can compute the final gas 
temperature, T(io, Aq). 

We use our numerical simulations to estimate 
the overdensity evolution, A(t,Ao). In Gadget the 
hydrodynamical equations are solved by the SPH 
method, hence, the fluid is represented by mass el- 
ements. The SPH-particles can be traced through 
the simulation. We group SPH-particles according 
to their overdensity at the end of the simulation, Aq, 
regardless if they are located in a galaxy, in a cosmo- 
logical filament or in the field. Then, we trace back 
the density of the particles and compute the average 
density of each group, see Fig.|4] For the integration 
of the thermal evolution of the IGM we approximate 
the overdensity evolution by an analytic expression, 
namely 



logA(t,Ao) = logAo 



G(a,Ao) 
G(l,Ao) 



(4) 



where we use the cosmic expansion factor, 
a = l/(z-|-l), as time variable. We approximate the 
time evolution by G(a, Aq) = |(a + l)™(i°gAo) _ i|^ 
and determine for each final overdensity the ex- 
ponent rn(logAo). We fit Eq. (|4]) to each curve in 
Fig.|4j where we take only the data for z < 3 into 
account. In the overdensity range — 2 < logAp < 6 
the resulting exponents are reasonably reproduced 
by m(logAo) = 1.3 log Aq - 0.26 (log Aq)^. This 
finally provides an analytic model for the overdensity 
evolution, A(t, Aq). 

In Fig. [5] the temperature-density relation, 
T{A,ta), obtained by integrating Eq. ([I]) using 
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Figure 5: Temperature-density relation, T(A,to) for the 
overdensity histories as given in Eq. Q (thick sohd hne). 
In addition two extreme scenarios for the overdensity his- 
tory are shown. First, all gas is compressed or decom- 
pressed to its final overdensity, Ao, at the very begin- 
ning, i. e. before reionisation takes place (dashed-dotted 
line). Secondly, all gas evolves at cosmic mean density 
until z — and is then instantaneously compressed or de- 
compressed to its final overdensity. The latter clearly re- 
sults in an adiabatic relation (dash-dotted line). Finally, 
the distribution of the SPH-particles in the simulation is 
shown. 



the analytic formula for the overdensity evolution 
described above, is shown. For comparison the 
density-temperature phase-space, as obtained from 
the numerical simulation, is also depicted. For high 
densities, A > 10"^, the analytic result, T{A,to), and 
the distribution of the SPH-particle coincides very 
well, because the thermal state is determined by 
thermal equilibrium. For low densities, most of the 
SPH-particles also follow T(A,to)- Finally, a signifi- 
cant fraction of the SPH-particles has temperatures 
substantially above T(A,to), since in the numerical 
simulation particles are also heated by dissipation at 
shock fronts. 

Which impact has the overdensity evolution? To 



answer this question we investigate two extreme 
scenarios. First, gas stays at mean cosmic density 
until z ~ and is then instantaneously compressed 
or decompressed to the final overdensity. This 
clearly leads to an adiabatic temperature-density 
relation, see Fig. [5] (dash-dotted line). Secondly, gas 
is brought to the final overdensity before reionisation 
takes place and has then constant overdensity. 
The resulting temperature-density (dashed line) 
relation is only slightly different from the result for 
the average cosmic density evolution (solid line); 
the transition to thermal balance occurs at lower 
overdensities, since the compressed gas has more 
time to cool. 

For the overdensity history introduced in Eq. ^ 
we find for low overdensities, A < 100, an effective 
equation-of-state, i. e. the temperature-density rela- 
tion follows a power-law, 



reos(A,t = <o) -ToA", 
with Tn = 3.85 x 10^ K and a = 0.59. 



(5) 



By plotting the cooling times in the temperature- 
density phase-space we can clearly see at which over- 
densities the gas is in thermal balance, see Fig.[6j For 
low overdensities, A < 10^, cooling has only a negli- 
gible effect, i. e. the thermal state of the gas is deter- 
mined by the cumulative photo- heating. In contrast, 
at high overdensities, A > 10^, the temperature lies 
exactly at the boundary between heating and cooling, 
i. e. it is in thermal balance. Therefore, in the range 
10^ < A < lO'^ the transition from the pure photo- 
heating to the equilibrium between heating and cool- 
ing takes place. 

4.3 The halo density profile 

We now set up the second part of our toy model, 
namely, the spherical halo model, with the aim to 
estimate the innermost gas density as a function of 
halo mass. 

For a spherically symmetric halo in hydrostatic 
equilibrium the pressure force and the gravitational 
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Figure 6: Theoretical result for r(A) with cooling (solid 
line) and without (dash-dotted line). The result without 
cooling follows the power-law T = 3.85 K {5 + 1)0-595 
Cooling times are indicated by colours and contour lines. 
Contours from left to right are at lOx Hubble time 
(dashed line), Hubble time (solid line), and O.lx Hubble 
time (dashed line). Cooling times are also indicated by 
colours, the blank area indicate the regime of UV heat- 
ing. The hashed area highlights the transition region from 
cumulative photo-heating to thermal balance. 

force are in balance, 

1 dP _ GM, 
p dr 

where G is the gravitational constant and Mtot(< r) 
is the total mass enclosed by the radius r. Assum- 
ing an effective equation-of-state, as given in Eq. (H), 
we can write the left-hand side of Eq. Q in terms of 
the gas overdensity. To this end we expand the gas 
density, p — An^mp, where rrip is the proton mass. 
The pressure can be written as P = TiiUikBT — 
(Anb//i) fcBToA", where we have used the molecular 
weight, /i = n\y/'E,ini, and Eq. ([s]). Now we rewrite 
the left-hand side of Eq. ^ , 



tot(< r) 



(6) 



p.mp 



(a + 1) A"- 



dA 
dr 



We simplify the right-hand side of Eq. ([6| by as- 
suming that the distribution of the total mass is de- 
termined by the dark matter. For the density profile 



of the latter we follow Ref. 45 



(8) 



with X = r/rg. The scaling radius, Tg, is related to 
the virial radius, rvir, of a halo by Ts — ry-n-jc, where 
c is the concentration parameter. The mass enclosed 
by the scaled radius x is 



M(< x) = l&TTrlp^ |ln(l 
= mTTrlp,F{x). 



(9) 



Evidently, the virial radius encloses the virial mass, 
Af^i, = M(< c). Dividing Eq. ^ by the virial mass 
leads to 

M(< x) = Mvir F{x)/F{c). (10) 

Using r — r^-^xjc and the expression Eq. ([7| we can 
rewrite Eq. (|6]) 



/imp 



(a + 1) A"-i dA 



GAUi,F{x) c 

Tvir F{c) X^ 



dx. (11) 



We can further simplify the expression on the right- 
hand side by introducing the maximum circular ve- 
locity. The circular velocity is given by 



IGM{< r) 



'GMvir c F{x) 



(12) 



X F{c) ' 



The quotient F{x)/x is independent from the halo 
properties and has a maximum value of ^ 0.22. 
Therefore, we can write for the maximum circular 
velocity 



GMvir 0.22 c 



(13) 



Eq. (Ill is now independent from the concentration 

F{x) 



(7) parameter of the halo, namely 



^(a + l)A"-idA=" 
/imp 



2 

max 

022 



dx. 



(14) 
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To determine the innermost gas overdensity we inte- 



grate Eq. ( 14 ) starting from the periphery of the halo, 
X — Xout, to the center, x = 0. Since the halo gas is 



highly ionised we adopt /i 
tral overdensity is 



0.59. The resulting cen- 



centre 



f To 



1) 



VlO^K 
ln(l 



\5.6 kms~i / 



1 - 



'^OUt ) 



where Aout is the gas overdensity at a;out- We 
estimate the outer overdensity by assuming that it 
is identical to the dark matter overdensity. 



From Eq. ([8]) we get 

Pdm 



A, 



(Pdi, 



(Pdm) X{1 + xy 



(16) 



The virial radius is determined by the cumulative 
overdensity, Ac{z), at which virialisation is expected. 
We use the approximation of Ref. [s] for Ac(z). For 
the virial mass we write now 



(17) 



where we have used Eq. ([9| for the last equality. Sub- 
stituting rs/{pdm) we can write for Eq. (16) 



Adm(a;out) 



Ac(z) 



F{c) a;out(l + Xout)^ ' 



(18) 



which allows us to estimate Aout in Eq. (15). 



We wish to give the central gas overdensity as a 
function of the maximum circular velocity. To utilise 
Eq. ( 18 1 we have to find a relation between Wmax and 



the concentration parameter. Maccio et al. [39] give 
an relation between c and the virial mass of a halo 
for z = 0, 



c(Mvir) = 10.47 X 



1012 h-^Mp 



-0.109 



(19) 



From our simulations we find that M^ir and the max- 
imum circular velocity are on average related by 



1012 h-^Mp 



^max 

210kms-i 



(20) 
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Figure 7: Central gas overdensity as a function of the 
halo circular velocity. We compute the central overden- 



sity according to Eq. ( 15 1 (thick solid line). Moreover, we 
plot the overdensity at the scaling radius, rs (thin solid 
line). The grey area indicates the transition region from 
pure UV heating to thermal balance, see Sec. |4.2| or Fig.|6] 
The blue dashed line indicates the characteristic circular 
velocity obtained from the numerical simulations. The 
thick and the thin red dotted lines indicate a condensa- 
tion criterion of A = 900 and 300, respectively. 



Moreover, the concentration parameter is exp ected to 
scale with redshift according to c ~ (1 -|-z)~i 



69 



73 



Thus, we obtain an relation between concentration 
parameter and Wmax, namely 



*^('^inax) 

10.47 



'^max 

210kms-i 



-0.32 I 



(21) 



which allows us finally to evaluate Eq. ( [Tsl ). 

Finally, we compute the central gas overdensity as 
follows: For a given circular velocity we determine 



the concentration parameter of the halo by Eq. (21) 



Then we estimate the outer gas density, here we 



assume that A(xout) ~ Adm(a;out), see Eq. (18), at 
three times the virial radius, Xout = 3c. Finally, we 
compute the central gas overdensity, Eq. ([15]), as a 
function of the circular velocity, see Fig. [7] 
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circular velocity at z = we choose /corr — 2.5. 



To estimate the characteristic circular velocity by 
our analytic model we have to give a criterion at 
which central overdensity gas condenses in a halo 
centre. The transition from the effective equation-of- 
state to the thermal balance takes place in the over- 
density range of 10^ < A < 10'^. For an overdensity 
of 900 our analytic model would match the charac- 
teristic circular velocity obtained from the numerical 
simulations. However, the resolution in the numerical 
simulation is not sufficient to reach the innermost gas 
density as assumed in our model. To mimic the nu- 
merical resolution we integrate Eq. ( |15[ ) only up to the 
scaling radius, see thin line in Fig.[7| In this case an 
overdensity of 300 corresponds to the characteristic 
circular velocity. This leads to two conclusions: First, 
the analytic model matches the numerical result very 
well. Secondly, better resolved simulations may re- 
sult in a lower characteristic velocity. However, gas 
mixing in the halo centre may generally prohibit an 
ideal hydrostatic equilibrium. 

4.4 Evolution of the characteristic 
mass 

In the previous section we have shown that our model 
well reproduces the characteristic circular velocity at 
z = 0. We can easily extend the model to other 
redshifts by determining the effective e.o.s and the 
criterion for condensation as a function of redshift. 
We derive the time-dependent e.o.s. by integrating 
Eq. ([T]), see Sec. 4.2 To obtain a time-dependent con- 
densation criterion we determine the position of the 
maximum in the density-temperature relation, see 
the solid line in Fig. [5] To compensate for the dis- 
crepancy between Amax and the actual overdensity 
for condensation we introduce a correction factor. 
To match the characteristic circular velocity obtained 
from the simulations this factor also comprises the ef- 



fect of the limited numerical resolution. Eq. ( 15 ) can 
be now rewritten 



V5.6kmS-V -^maxJcorr 



To a + 1 



104 K a 



(22) 



where we have neglected the contribution of the 
outer halo boundary. To match the characteristic 



Eq. ( 22 1 indicates that the characteristic circular 



velocity is essentially determined by the density- 
temperature relation at any time. All parameters 
on the r.h.s., namely slope and normalisation of 
the e.o.s. and Amax, are derived from T(A,t). 
Fig.js] (upper row) shows how the characteristic 
circular velocity is composed of three parameters 
of the temperature-density relation. The upper left 
panel shows T(A) for several redshifts. The next 
three panels show the evolution of the parameters 
of T(A,i): The temperature To, the overdensity 
at the maximal temperature, Amax, and the slope 
dlogT(A,t)/dlogA at A = 1. The right panel 
shows the resulting Ui/2: which nicely matches the 
characteristic circular velocity obtained from the 
numerical simulations. In the range 3 > z > the 
circular velocity is virtually constant. This is the 
result of the decreasing Tq and the increasing Amax- 

The redshift evolution of the characteristic mass 
can be derived from characteristic circular and 
Eq. (131 and Eq. (17). Combining these equations 
and using (p) (x {z + 1)^ we find 



(23) 



This is in agreement with the results found earlier 
More precisely, the dominating terms in M1/2 



22 



are introduced by the redshift dependence in the def- 



inition of the virial radius, Eq. (17 1 



4.5 Variations in the UV background 

As shown above the cumulative heating of the IGM 
by the UV background prevents the gas cooling in 
the centre of small mass haloes. It is therefore evi- 
dent that the characteristic mass depends on the UV 
background and its evolution over cosmic times. We 
use our model worked out above to assess the im- 
pact of modifications of the UV background model 
on the characteristic circular velocity. We first con- 
sider changes of the UV flux density and secondly of 
the UV spectrum. 

UV flux density The flux density of the UV back- 
ground has only a small effect on the temperature of 
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Figure 8: Left panels: The density-temperature T(A, t) for several redshifts and three different UV-background 
models. The next three columns show the redshift evolution of the parameters of T{A,t): The normalisation and 
the effective e.o.s., To, the position of the maximum, Amax, and the slope of the e.o.s., a. Finally, for each UV- 
background model the resulting characteristic circular velocity is shown, see Eq. ( |22[ ). The upper right panel shows 
also the characteristic velocity obtained from the numerical simulations. 



the IGM. The reason is that a higher flux reduces 
the fraction of neutral hydrogen in the tenuous IGM, 
hence, the higher flux flnds fewer neutral atoms to 
ionise. As a result the heating rate is essentially inde- 
pendent from the flux density, when photo-ionisation 
dominates over coUisional ionisation and as long as 
the spectrum remains unchanged. In this case the 
fraction of neutral hydrogen, xhi = 1 — a^Hii, can be 
derived from the balance between photo-ionisation 
and recombination 



(24) 



The IGM is generally highly ionised, hence xhii ^ 1- 
In this case the electron density is proportional to the 
hydrogen density, rig ~ nn- As a result the heating 
due to hydrogen ionisation becomes 



chii 
i 7HI 



(25) 



The heating therefore depends only on the ratio 
eHii/r7Hi, i- e. only on the average energy input per 
ionisation event. In the limits made above the heat 
input is independent from the actual flux density. In 
accordance with this result we have found in Ref. 22 



that the characteristic mass depends only httle on 
the UV background flux density. Even very strong 
variations of the flux causes only a small shift of the 
characteristic mass, showing that the latter is sensi- 
tive only to the thermal state of the IGM which is 
httle affected by the actual flux density. 

UV spectrum Significantly more important for 
the thermal evolution of the IGM is the spectrum of 
the UV background. As a first guess one may expect 
that the final temperature of the IGM is directly 
proportional to the energy per ionising photon. 
However, since the recombination rate, aHii(2^), 
decreases with increasing temperature the effect is 
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overdensity A = p /(p) 

Figure 9: Theoretical results for temperature distribution in the IGM as a result of different UVB heating models. 
Left panel: The thick solid line indicate the result for our standard model according to [20| . Color coded is the 
cooling time. The blank area indicates the region where the gas is heated. Three contour lines are shown, they 
indicate from the left to right a cooling time lOx, Ix, and O.lx the Hubble time. The dashed blue line indicates a four 
times increased energy per ionising photon only for z > 2. Middle panel: Results for four times increased energy per 
ionising photon for all redshifts. For comparison the dashed blue line indicates temperature distribution resulting 
from the standard model as shown in the left panel. Right panel: Results for ten times increased energy per ionising 
photon for all redshifts. For comparison also the eight times increased energy and the standard heating as shown in 
the left panel are indicated. 



not that strong. Fig. |9] shows the results for increased 
energies per ionising photon compared to the 20 
model. We have computed models with four times 
and with eight times increased energy per ionisation 
event. The higher average energy shifts the effective 
equation-of-state to higher temperatures. Also the 
shape of the temperature-density relation changes. 
The maximum in the relation shifts to higher 
densities. In particular, the maximum changes 
significantly between the four times and eight times 
increased model. The effective equation-of-state runs 
in the cooling regime virtually parallel to lines of 
constant cooling times. When Tq > 1.5 x 10^ K these 
cooling times are above the Hubble time. In those 
cases the cases the transition to thermal balance 
appears at very high ovcrdcnsities. 

We use the models introduced above to estimate 
the impact of the UV spectrum on the characteristic 
circular velocity. To this end we compute the evolu- 
tion of the characteristics of the temperature-density 
relation, see Fig.|8] For the four times increased en- 
ergy per ionising photon we find that the characteris- 



tic circular velocity amounts to ~ 60km s~^ for z <2. 
For the eight times increased model the characteristic 
velocity get very high, since Amax shifts to very high 
overdensities. 



5 Summary 

We discuss high-resolution simulations including 
radiative cooling and photo-heating. We focus on 
the baryon content of dwarf galaxy-sized dark matter 
haloes. To this end, we have simulated the structure 
formation in two distinct cosmological environments 
and we found that the characteristic mass below 
which haloes are baryon deficient does not depend on 
the environment. Moreover, we found that in haloes 
in the mass range of ~ 10^ to 10^° Mq in general 
gas condensation stops at some redshift between 
reionisation and today. In our simulation those 
haloes are able to retain the already condensed gas, 
allowing further star formation. We have derived the 
characteristic mass in newly accreted matter, which 
is somewhat higher than the characteristic mass of 
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the baryon fraction in the haloes. The former is 
more appropriate for semi-analytic modelling. 

Tinker and Conroy found by a Halo Occu- 
pation Distribution analysis that the magnitude of 
galaxies which reside in haloes with a mass about 
the characteristic mass is Mr ^ —10. Moreover, 
they found that more luminous void galaxies, i. e. 
Mr > —10, are preferably located close to the 
walls of the voids. Therefore, photo-heating would 
reduce the baryon content of the small-mass haloes 
populating the entire void volume, rendering these 
galaxies even more dark. Current surveys which 
cover large void volumes are not sensitive enough to 
detect the galaxies affected by photo-heating. 

The failure of condensing gas can be traced back 
to the cooling time in the halo centre, which may 
exceed the Hubble time. We have set up a spherical 
model for the haloes gas density profile. Crucial for 
this model is the effective equation of state which we 
have derived by integrating the thermal history of 
the IGM using an approximation for the overdensity 
evolution as a function of the final overdensity We 
have justified this model by considering the density 
contrast evolution in the simulation. This spherical 
model allows us to derive the characteristic mass for 
different heating histories of the IGM. 

We find that the mass at which galaxy formation 
seems to fade can be explained by a heating model 
which incorporates 6-8 times more energy per 
ionising photon than given in the Haardt & Madau 
model 
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For redshift z > 3 it is known from 
analyses of the Lymana forest that the temperature 
of the IGM is higher than obtained by a naive 
application of the Haardt & Madau model. The 
derived temperature implies a 3-4 times increased 
energy per ionising photon, or a more efhcient 
heating due to non-equilibrium effects or radiative 
transfer effects. We found that if the average 
temperature of the IGM is today as high as > 10"* K 
dwarf galaxy formation would be suppressed at a 
mass scale consistent with that derived, e.g., by the 
conditional luminosity function. Therefore, the mass 
scale for suppression of dwarf galaxy formation may 



be considered as a measure for the temperature of 
gas in the surroundings of dwarf galaxies. The heat 
source could be the UV background or alternatively 
any feedback of the galaxy. 
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